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Abstract — To account for the randomness of propagation 
channels and interference levels in hierarchical spectrum sharing, 
a novel approach to multihop routing is introduced for cognitive 
random access networks, whereby packets are randomly routed 
according to outage probabilities. Leveraging channel and in- 
terference level statistics, the resultant cross-layer optimization 
framework provides optimal routes, transmission probabilities, 
and transmit-powers, thus enabling cognizant adaptation of 
routing, medium access, and physical layer parameters to the 
propagation environment. The associated optimization problem 
is non-convex, and hence hard to solve in general. Neverthe- 
less, a successive convex approximation approach is adopted 
to efficiently find a Karush-Kuhn-Tucker solution. Augmented 
Lagrangian and primal decomposition methods are employed 
to develop a distributed algorithm, which also lends itself to 
online implementation. Enticingly, the fresh look advocated 
here permeates benefits also to conventional multihop wireless 
networks in the presence of channel uncertainty. 

Index Terms — Routing, cross-layer optimization, multihop 
wireless networks, cognitive radios, random access, channel 
uncertainty, convex approximation, distributed computation. 

I. Introduction 

Resembling traditional routing protocols for wired net- 
works, their counterparts for wireless networking generally 
utilize optimization tools such as shortest path routing to 
find optimal route(s) based on the network connectivity graph 
abstraction fTJ. Early on, links among nodes were quantified 
based on a disk model capturing only distance-based determin- 
istic losses. Upon recognizing the inadequacy of disk models 
for the broadcast wireless interface [2|], a weighted graph 
accommodating more sophisticated performance metrics was 
adopted; see e.g., 0, 0, and the stochastic routing approach 
in 0, where link weights capturing packet delivery probabili- 
ties were exploited to develop optimal routing schemes. These 
schemes are particularly attractive for energy-limited nodes, 
primarily because the resulting routing strategies promote links 
with higher reliability, thus decreasing the number of packet 
lost due to fading 0. 

In a hierarchical access setting, interference levels can 
not be acquired accurately due to the lack of explicit inter- 
system cooperation |7j. As a result, random shadowing and 
small-scale fading effects, along with dynamically changing 
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activities of licensed users, accentuate the uncertain nature of 
wireless cognitive radio (CR) links. The effects of random 
interference on CR links from primary user (PU) transmitters 
is called upon in 0, where source-to-destination paths that 
are most likely to meet prescribed end-to-end requirements 
are found based on predicted link capacities. Leveraging the 
situational-awareness provided by spectrum occupancy detec- 
tion schemes, a graph whose link weights reflect the amount 
of spectral resources available per CR-to-CR link is employed 
in 0, where optimal routes are obtained via Dijsktra or Bell- 
man Ford-like algorithms. A two-phase approach is proposed 
in iflOl . where nodes in the network first obtain an expected 
route cost and a set of candidate forwarding nodes, and then 
route traffic across paths with higher spectrum availability. 
In irTT1l . the average link availability is invoked to develop 
a routing scheme that avoids network zones with unstable 
CR connectivity. Link availability in ifTTI is computed in a 
probabilistic sense based only on the statistics of primary 
user (PU) activities. A PU coverage map supplied by sensing 
schemes is employed in 1121 to identify spectrum opportunities 
in space, and devise routing strategies supporting multiple 
classes of CR quality-of-service (QoS) demands. 

The aforementioned works offer valuable insights on route 
formation and management based on the average availability 
of CR links, and predicted link capacities. However, in a 
hierarchical access setup, link capacities are unknown and may 
change abruptly because of time-varying PU activity patterns, 
dynamic shadowing, and diverse QoS constraints. In this 
context, a cross-layer design approach to obtain both optimal 
routes and physical and medium access parameters that dictate 
the packet forwarding capabilities is therefore well motivated. 
To this end, the present paper exploits propagation channel 
statistics to develop a statistical routing approach whereby 
nodes not only compute optimal routes, but also optimal 
link reliabilities by controlling transmit-powers and medium 
access control (MAC) parameters. The novel approach ac- 
counts explicitly for the randomness of propagation and the 
medium access interface, to allow spectrum-cognizant routing 
of data packets, while enforcing PU interference protection 
(Section Hill. 

In spite of the non-convexity of the associated cross-layer 
optimization problem, a successive convex approximation is 
pursued to find a Karush-Kuhn-Tucker (KKT) solution effi- 
ciently (Section [TTTb - Enticingly, feasibility guarantees offered 
by the successive convex approximation algorithm naturally 
suggest an online implementation of the algorithm whereby 
nodes do not necessarily wait for the successive convex 
approximation iterations to converge, but rather use network 
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parameters as they become available. 

However, the communication overhead incurred to acquire 
channel statistics at a central node, and subsequently dissem- 
inate optimal network parameters can become prohibitive as 
the network size increases. To alleviate such a message-passing 
burden, and address scalability and robustness concerns, a dis- 
tributed algorithm is also developed by invoking the alternating 
direction method of multipliers and the primal decomposition 
method (Section |IV| ). Finally, suitable conditions are estab- 
lished to ensure that packets are eventually delivered to their 
destination when routes, medium access and physical layer 
parameters are regularly updated to track channel statistics 
and topology dynamics (Section [Vb. 



A. Preliminaries and problem formulation 



N 



Consider a wireless CR network with N nodes {U n } 
sharing spectral resources with an incumbent PU system 0. 
Leveraging the spectrum awareness provided by spatio- 
temporal sensing schemes (13l , fBl . CRs collaborate in 
routing data packets to a sink node E/jy+i, while respecting the 
PU-CR hierarchy. The CR network is modeled as a digraph 
to account for the possible lack of link bi-directionality. The 
dynamic and stochastic nature of the CR propagation ambi- 
ence, along with the possibly minimal amount of topological 
information motivate consideration of random medium access, 
as well as stochastic routing strategies J4), J5), lfl5l . In this 
context, a CR node U n transmits with probability /i„ £ [0, 1], 
and decides whether to route packets toward a neighboring 
node Ui with probability t n ^i £ [0,1] per time slot. As 
packets are forwarded to neighboring nodes according to 
probability mass functions, it holds that J2i^ n tn-yi = 1> f° r 
alln= 1,...,N. 

Communication of data packets over a wireless network 
depends not only on transmission and forwarding decisions, 
but also on the intended link reliability. In case of unsuccessful 
packet decoding due to fading- or interference-induced link 
outages 121, a packet not eventually routed by U n will remain 
in Un & queue, and its transmission will be re-attempted in 
a subsequent time slot (possibly to a different neighboring 
CR). To capture channel- and interference-induced sources of 
uncertainty, let r n _>.j £ (0, 1] denote the probability that a 
packet transmitted from node U n is correctly decoded (and 
thus successfully received) by Ui. 

Assuming that link reliabilities {r„_^} are known by, e.g., 
computing the packed error rate of preceding sessions, a 
stochastic routing framework for maximizing users' exoge- 
nous rates was introduced in 15] . However, because of the 
volatile CR channel characteristics, time-varying PU activity 
patterns, and diverse QoS constraints, {r„_j.i} may change 
abruptly during the network operation. Hence, {r n ^} may 
not be known in advance. Building on first- and second-order 
statistics of the PU interference, as well as those of node-to- 
node channels, a statistical routing approach yielding optimal 
(i) routes, (ii) transmission probabilities, and (iii) transmit- 
powers is put forward in the ensuing section. 



II. Statistical Routing Framework 

Data percolation through a wireless network is captured by 
the product packet delivery probabilities r„^i}. When 

random access is employed as MAC, it is common to consider 
a packet lost when collisions among CR transmissions occur. 
With I n i denoting the set of nodes whose transmissions 
interfere with link U n — > Ui, the probability of collision-free 
packet transmission from U n to Ui is given by Iljei ■ (1 — Mj)- 
A widely-accepted criterion for successful packet reception 
is to require the signal-to-interference -plus-noise ratio (SINR) 
to stay above a certain threshold |f2), J6], which is generally 
determined by the receiver structure, modulation, and coding 
scheme. Let g n ^i denote the channel gain between U n and 
Ui, modeling the effects of path loss, log-normal shadowing, 
and Nakagami-m small-scale fading lfl6l . Then, the SINR of 
link U n — > Ui can be expressed as 



In 



„1 i Y^-^s ^ 



(1) 



where af stands for the receiver noise power at Ui, p n £ 
(0,p" lax ] denotes the transmission power of U n ; and irs.i the 
interference perceived from PU transmitter S = 1, . . . , Ng. 
Randomness of {j n ^i} in (H) emerges due to the shadowing 
and small-scale effects on the PU interference {tts^}- Fur- 
thermore, CR-to-CR gains {g n ^i} may be known imperfectly 
because of insufficient time for channel training. Nonetheless, 
CR-to-CR and PU-to-CR deterministic path losses, and statis- 
tics of shadowing and small-scale fading can be acquired and 
used. To this end, it is useful to recall that the distribution of 
channel gains {g n ^i} can be approximated as log-normal lfl6l 
Ch. 2], ifTTll . Furthermore, the Fenton- Wilkinson result jl8| 
asserts that the distribution of SINRs {7n-n} in (dJ can be 
well-approximated as log-normal too, with mean and variance 
expressed in terms of the first- and second-order moments 
of {gn^i} an d {i"s,i}; see IfTTl for a detailed derivation. 
Consequently, r n _>j := 10 log 10 7 n _>j will be approximately 
Gaussian distributed with mean P n + m n ^i, where P n := 
101og 10 p„, and variance denoted by cr^^t- The probability 
r„_s.i that a packet transmitted from U n is correctly received 
by Ui can thus be expressed as 

r n ^i = [ (1 - /-ij)Pr{7„^ l > 7„^ l } 

« n i 1 - q - p "~ m "^ (2) 



where Q(x) := f°° —?=e ^ dx is the standard Gaussian 

^ V 1 _ Jx V27T 

tail function, j n ^i is a prescribed SINR threshold, and 
r„-H := 10 log 10 TW-H- Similar to fl9l . the main interest 
here is in the tail of the complementary cumulative density 
function (ccdf) of the SINR; in this case, the Fenton- Wilkinson 
method is known to provide accurate approximations for all 
the propagation scenarios of practical interest lfT6l Ch. 3], ATI . 
11201. 

Using (O, the link reliabilities {r n ^i} can be expressed 
in terms of the MAC variables {^„} and the physical layer 
quantities {P n , m„_H, cr^J. Therefore, with {m„^ i5 <r n ^ t } 
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Fig. 1. Input and output flows at node C/ n under queue stability. 



known parameters, the optimal routing strategy will be ob- 
tained by optimizing over {p n , t n ^i, Pm}- The next step is to 
model exogenous data packet arrivals at U n from its applica- 
tion layer by a stationary stochastic process with average rate 
p n G (0,1] per time slot. Suppose also that each CR node 
maintains a backlog to cache exogenous and endogenous^ 
packets that have to be routed toward the destination Un+i- 
Aggregate queue service rates depend on the joint queue 
occupancy distribution. This results in a generally asymmet- 
ric system of interacting queues, whose stability region is 
challenging to analyze even for simple systems. Nevertheless, 
assuming as usual fully backlogged queues per node (21] 
yields a sufficient condition for queue stability that can be 
conveniently used as a constraint in rate-oriented routing 
optimization. In the advocated dominant system, users with 
empty queues transmit "dummy" packets, and consequently 
queue sizes are never smaller than those in the original system, 
if both systems start from the same initial condition. 

Let A„ denote the average aggregate rate of endogenous 
packet arrivals at U n , which coincides with the rate of packet 
departures if queues are stable. Then, queue stability implies 
that {p n }n=i an d {^n}n=i abide by the flow conservation 
constraints 0, E2 (cf. Fig. [Tj 



>j ^ ^iti-+'nXi-+nt V (3) 



where A/" n _> := {j\r n ^j > 0,j = 1,...,N + l,j ^ n} is 
the set of nodes that decode J7„'s transmissions with non-zero 
probability, and A/"_s.„ := {i|r.;^„ > 0, i = 1, . . . , N, i ^ n} 
the set of nodes that route packets through U n . For queue 
stability, Loynes' Theorem |23l asserts that for stationary 
arrival and departure processes (the latter are stationary in 
the dominant system) a sufficient condition for stability is 
A n < Pn, for each CR U n ; and a necessary condition for 
stability is A„ < /i„ (cf. Fig. [TJ. 

To complete the formulation, consider Nr actual or po- 
tential PU receivers, whose locations have been estimated 
via sensing l24l . and let tg ax denote the maximum average 
interference that can be tolerated by PU receiver R Q, 1251 . 
Further, let Mr C {U n }„ =1 be the (sub-)set of CR nodes 

'"Exogenous" packets of a CR node are those generated from its application 
layer. On the other hand, "endogenous" packets refer to those received from 
the neighboring nodes of a CR node, and are to be routed by the network 
layer; see also Fig.fT] 



located in the proximity of PU R (not necessarily the entire 
CR network, as some CRs may be sufficiently far apart and do 
not interfere with PU R). Transmissions by CR U n undergo 
random shadowing and small-scale fading effects before arriv- 
ing at close-by PU nodes. Approximate the channel gain g n ^R 
between CR U n and the PU R as log-normal distributed |[T7l . 
and define a binary random variable a n G {0, 1}, independent 
of g n ^R, taking the value 1 with probability p, n , and with 
probability 1 — /i„. Then, the average interference experienced 
at PU R is given by (k := 0.1 ln(10)) 



L R := E < a n p n g n ->R 



N 

,, „K-P7i+Kmn^H+V CT „_+H <T .max 
pnC- ~ ^ L R, ■ 

n=l 



(4) 



Variables {P n , p n , p n , A„}, and {t n ^i} satisfying the con- 
straints (f3]l and (0]i can be supported by the wireless CR 
network. It is certainly desirable to design the network by 
selecting a feasible set of variables that are optimal in some 
sense. To this end, consider a concave utility U n (p n ), and a 
convex cost C n (P n ), representing the reward of exogenous rate 
p n and the cost of power P n for node U n , respectively 1261 . 
Notice that p n is the average rate of packets generated at the 
application layer of node U n to be eventually delivered to 
the sink Un+i l22l . (26); thus, p n represents an end-to-end 
performance metric. Capitalizing on the statistical description 
of SINRs and CR-to-PU channels, the statistical routing 
problem is formulated as: 



N 



N 



(PI) max Y U n (pn)~ VC„(F„) (5a) 



{/i„>0},{A„>0} 
{t„^i>0} 

subject to 



Pn+ 




n-yj 






Vn = 


1,.. 


.,N 


(5b) 


E 


t n -+i < 1, Vn = 


1,.. 


.,N 


(5c) 




A„ < p n , Pn < 1, Vn = 


1,.. 


.,N 


(5d) 




P n < P™. Vn = 


1,.. 


.,N 


(5e) 




max \ \ KP„+sm„ 
L R -± ^ Pn e 

n<£Af R 


-^R + 


3 






VR = 


= 1,- 


■-,Nr 


(5f) 



with {r n ^j} given by ©, and P™ ax := 101og 10 p™ ax . 

The non-convexity of constraints ( |5bl and ( fSfb makes prob- 
lem (PI) non-convex, and thus hard to solve. Furthermore, 
function Q( ) in (O is difficult to handle in an optimization 
problem. In the next section, an approximate but efficiently 
solvable version of (PI) will be formulated. But first, some 
remarks are in order. 

Remark 1 (Monotonically non-decreasing utilities). It follows 
from ll27l Thm. 5], that (PI) is optimally solved by setting 
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{A„ = fi n } if each utility function U n (p n ) is monotonically 
non-decreasing. As many practical utilities satisfy this condi- 
tion, U n {pn) will be hereafter assumed non-decreasing, and 
variables {A„} will be dropped. Strictly speaking, the choice 
{A„ = fx n } will lead to a solution of (PI) where queues 
may or may not be stable l23l . On the other hand, condition 
A ra < p,„ is challenging because it entails an open constraint 
set. From a practical perspective, queue stability can be readily 
ensured by imposing in d5eb the condition A„ + e < p n , with 
< e <C 1 small enough, and replacing variables {p n } with 
{fin := fi n - e}. □ 
Remark 2 (Conventional multi-hop networks). The proposed 
routing framework can be considered also for non-CR mul- 
tihop random access networks when node-to-node channels 
can not be estimated accurately - what could emerge with 
e.g., a mobile ad hoc topology. Optimal routes and link 
reliabilities can be obtained by solving (PI), after discarding 
the interference constraints OB . and re-defining the signal-to- 
noise ratio (SNR) of link U n -> Ui as j n ^i = p n g n ^i/(rf, 

□ 

Remark 3 (MAC protocol). Since a random access protocol 
is adopted, a packet is as usual deemed lost when collisions 
among CR transmissions occur [cf. ©], and no mutual in- 
terference is explicitly modeled in (|T). However, the solution 
approach presented in the ensuing section can be effectively 
employed when different MAC strategies such as, e.g., car- 
rier sensing medium access and orthogonal transmissions are 
utilized by the CR nodes; see also ll26l . □ 

III. Tractable Routing Protocol 

To convexify constraint Oft it suffices to consider the loga- 
rithmic change of variables fl n := ln(/i n ). As for the flow con- 
straint (l5bl . consider first introducing auxiliary variables {f n } 
representing the probability of CRs to remain silent, together 
with the extra constraints \i n + v n = 1, for n = 1,...,N. 
Further, a simple way to obtain a tractable approximation of 
Q(x) consists in exploiting the commonly used upper and 
lower bounds proposed in J28), 11291 . which are very tight for 
x > V2/2. Taking advantage of these bounds, and performing 
again a logarithmic change of variables v n := ln(^„), the 
probability can be (tightly) bounded as 

> ** x (1 - _L e -K Pn+m C-: I r ""') 



4 



< e ^ 3 ez,„ v i x ( 1 - a x 



(6) 



(7) 



where a\ = 0.28, and = 0.64 |29|. The premise for 
adopting the aforesaid bounds is that the decoding rate of CR 
links is at least w 0.7. This condition is met in practice if CRs 
and PUs are sufficiently far apart (see, e.g. Ifl9l ). Furthermore, 
maximum packet error rates required for data and speech 
transmissions are generally considerably lower than 0.3 J6). 

Consider now using the upper bound (0 for the incoming 
traffic, and the lower bound (|6]i for the outgoing flows. As (|6]l 
and © are tight, this replacement not only yields a tractable 



optimization problem, but also does not sacrifice optimality of 
the outcoming rates. With the logarithmic change of variable 
i„_M = ln(£ rl _>i), and after introducing auxiliary variables 
{iln-¥i > and {y n ^i > \/2/2}, constraint ( Bbl can be 

approximated as 



p n + £ /<" +t -+^^ 



12 



1 



E ■ 



< 



+ E « 

with the auxiliary constraints 

1 

a^iO/n^i) < P n + m n ^i - T n ^i (9) 

vLniVi^n) >Pi + mi^n - f (10) 

For notational convenience, define the variable vector x„ :— 
{y 3 ^ ni y n ^i}] T per node U n , for n = 
1,...,N. Upon re-expressing in a compact form the flow 
constraint ([8]) as /„({x„}) < 0, and defining the constraint 
set B Xn per node U n as 



Pn < P" ax ,and (11) 



where (non)negativity of the variables is left implicit, problem 
(PI) can be re-formulated as 



N 



N 



(P2) max 



J2U n ( Pn )-J2 C n(Pn) (12 

^ 1 „-, 1 



a) 



subject to /„({x„}) < 0, n = 1, 

E] /'n,™( x «) — t i? aX ' -R = 1) 



,N 
,N R 



(12b) 
(12c) 



n<£j\f R 



where / tR ,„(x„) 

Constraints ([8]l are still non-convex because the last two 
sums (with their signs) are concave, and likewise © is also 
concave. Nevertheless, the structure of (P2) allows convex 
approximation methods for obtaining its solution efficiently. 
Among candidate methods, the successive convex approxi- 
mation approach ll30l is well suited for the problem at hand 
because it guarantees first-order KKT optimality under mild 
regularity conditions. 

A. KKT solution via successive convex approximation 

The general successive convex approximation method is 
outlined first. Suppose that the objective function to be max- 
imized is concave in the optimization variables x, and the 
constraint set is the intersection of a set A := {x|/ n (x) < 
0, n = 1,...,N} with a convex set B, which captures 
convex constraints, if any. Assume that /„(x), n = 1, . . . , N, 
are differentiable but generally non-convex functions. Then, 
starting from a feasible point x^ ' e AnB, a series I = 1, . . ., 
of surrogate problems is solved, where A is substituted per 
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iteration i by a convex set A^'. Since the intersection of 
convex sets yields a convex set, the resulting optimization 
problems are convex. For each n = 1, . . . , N, let / n (x; x^) 
denote the surrogate convex function for /„(x), which may 
depend on the solution xW to the problem of the previous 
(I — l)-st iteration. Then, the convex set A^' is constructed 
as A w := {x|/„(x;xW) < 0,n = 1,...,N}. Provided 
that each function /„(x;x^'), n = 1,...,N, is convex, 
differentiable, and satisfies conditions 1301 

cl) /„(x)</„(x;xW), Vxe#nB 
c2) / n (xW)=/ n (xW;xW),and 
c3) V/„(xM)=V/„(x»;x«) 

the series of solutions to the approximate problems converge 
to the KKT point of (P2) ||30l . 

In order to apply the successive convex approximation 
method to (P2), surrogate constraints for the non-convex 
constraints must be determined. The first three terms in ^ are 
convex, whereas the fourth and fifth terms are concave. Letting 
_ e x!+px 2 -ax 3 re p resen t one f the non-convex summands, a 
convex surrogate function satisfying cl)-c3) can be obtained 
by replacing the non-convex summands with the affine func- 
tion 



a x 1 +fjX2 — ax 3 ^ ^x\ '+px) 2 ' — ax' 3 



Xl 



-/3(x^ - x 2 ) - a( X W - x 3 ) - 1 



(13) 



As for ©, an upper-bound of y/ ij n ^j can be obtained via 
the supporting hyperplane, and the resulting surrogate convex 
constraints become 



Un— >j 



Pr, 



f„^<0. (14) 



Overall, the problem to solve in the £-th iteration is given by 
(P2) with (O replaced by ( TT4T > to form the surrogate constraint 
set B n , and by employing ( [Oi l along the feasible points 
{xi;P}n = i to obtain a surrogate convex flow conservation 
constraint fn\{x n }) < 0; that is, 



N N 

max y~] U n {p n ) - y~] C n (P n ) (15a) 

J n—1 n— 1 



(P2^) 

subject to (HSDl and /^({x,,}) < 0, V 



(15b) 



Problem (P2 W ) is convex, and thus efficiently solvable using 
interior-point methods PP . It is worth mentioning that the 
solution of (P2^)), £=1,2,... always lies inside the feasibil- 
ity region of the original non-convex problem (P2) [30]. This 
observation suggests readily an online practical implementa- 
tion of the algorithm whereby node U n does not necessarily 

wait for the successive convex approximation algorithm to 

(£) 

converge, but rather relies on x„ as and when it becomes 
available. In the limit (i.e., for I > 1), x^ will be KKT- 
optimal. An online implementation of the iterative optimiza- 
tion allows tracking of slow variations in the network topology 
and SINR statistics. 



IV. Distributed statistical routing 

To obviate the high communication cost associated with the 
collection of channel statistics for all links at a central process- 
ing unit, and the subsequent dissemination of the optimized 
variables, it is of prime interest to solve (P2) in a distributed 
manner. A distributed cross-layer optimization algorithm is 
also desirable because of its scalability with regards to power 
requirements and network size, and robustness to isolated 
points of failure. 

Distributing (P2) is tantamount to developing a distributed 
solver for each of the convex problems (P2< £ )), I = 1, 2, . . .. 
To this end, it is necessary to decompose (P2^)) into smaller 
sub-problems, which can be locally solved by nodes {U n } 
via local message exchanges. Unfortunately, the interference 
constraints (1 1 2ct challenge decomposability, as they couple 
portions of the CR network. Furthermore, for each U n , 
constraint (1 1 5bb involves variables pertaining to the one-hop 
neighboring nodes Ui £ A/l^. n , and to CRs in the collision- 
related sets {Ijn} and {I m }. To overcome the first hurdle, 
consider first the following problem 



(P3 W ) 5({CW1):= 



N N 

max ^ U n (p n ) - ^ C n (P n ) 

subject to /„({x„}) < 0, Vn 

A«,n(x„) < Vn,R 



(16a) 

(16b) 
(16c) 



where the interference t'^ lax for PU R is pre-partitioned in 
given per-CR fractions {{^(t)}} n eM R - Problem (P2 W ) 
will be revisited later on. Then, collect local copies of 
Xj^ rl := [tj^n, Pj 7 jlj] T at node U n , for each j £ J\f-y„ 
into a vector Xj^n^ := [tj_j. nj „, Pj. n , P^j.n] T ■ Likewise, let 
{v n ,m} denote local copies of {v m \m £ l n }, with I n := 
(UieA^„^2ni)U(Uj e AT^ ?l X,' n ); i.e., local copies of v m for users 
that may interfere with [/„'s transmissions. Then, (P3^) can 
be equivalently re-formulated as 



N 



( P4») 



max 



N 

£ 

71=1 



C n {P n ) (17a) 



subject to 

/nfxnjx^},^})^ Vn (17b) 

/ lR) „(x n ) < t^(£), Vn,R (17c) 

Xj_> n = x.j^. ntn , j£A/^„,Vri, (17d) 

v m = iw, m£l n ,\fn (17e) 

where the notation f„(x n , {xj^.„. n }, {v m ,n}) emphasizes the 
dependence of the surrogate flow conservation constraint /„(•) 
on the newly introduced local variables. Problem (P4^) is 
amenable to a distributed solution, where dl7dt - dl7eb can be 
enforced by means of local message passing. 

Suppose that there is a non-zero probability (possibly multi- 
hop) directed path connecting U n to nodes U m £ Af-> n U 
Af n ^t U In', i-e-, nodes coupled in the optimization problem. 
If not, a control channel can be employed as usual. Problem 
may be solved in a distributed manner using the dual 
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Algorithm 1 Distributed algorithm for (P4^) 

Assumption: bidirectional links, or bidirectional control channels. 
Use solution of (P^ -1 ') to initialize variables, 
for 1 = 0,1,... (repeat until convergence) do 

Receive multipliers {q. n ,i(l)}iEM n ->, {^n,p(l)} £ e{p\u rl ei p } ■ 
Update x„ := {x„, {xj->„,„}, {v m ,n}} via dl9a| > 
Transmit x„_>i(7) to Ui G A/k-y, and x ■,_►„,„(£) to Uj G N^ n 
Transmit {v m ,n} to U m G X n via neighboring nodes. Forward 

{Vm,j}j£Af^ n tO Ui G Afn^- 

Receive Xj_> n (Z) from f/i G A/" n ->, and x n _>; from (7; G 
M,- 

Receive ^^(i) from Up€ {U p \U n G X p } 
Dual update via J19ct -( fT9dl . 
Transmit multipliers q n> j(i + 1) to Ui G N n ^>- 
Transmit v m ,„(l + 1) to £/ m G In via neighboring nodes, 
forward {v m ,j}jeM-, n + 1) t0 ^« e .Mi->. 
Use parameters x n to transmit data in case of on-line imple- 
mentation, 
end for 



sub-gradient method 1321 . 031 . However, recovering the pri- 
mal variables {x„} from the Lagrange multipliers optimizing 
the dual function is not always guaranteed if the objective 
in (I17ab is not strictly convex, and the step-size in the sub- 
gradient ascent is constant. Furthermore, primal averaging can 
not be performed in this case, unless the equality constraints 
are appropriately relaxed J34). 

One effective remedy is offered by the alternating direction 
method of multipliers (ADMoM), where the optimization 
argument in (P4^) is augmented with a quadratic regulariza- 
tion term corresponding to the squared norm of the equality 
constraints l35l Sec. 3.4]. Specifically, letting {q^n} and 
{v m ,n} denote the multipliers associated with the equality con- 
straints (I17dl i and (I17eb . respectively, the partial quadratically- 
augmented Lagrangian function is given by 

N 

£({x n },{q,-„},{n m ,„}) := ^2 \-Un(Pn) +C n (P n ) 



i >n,n || 2 



(18) 



where x„ := {x„, {Xj-y ni „}, {v m ,n}}, and c > is an 
arbitrary constant. Notice that £(■) is defined over the primal 
feasible region A := r\^ = iA n , with A„ '■= {x„|x„ G 
B Xr , (I17bl i. d!7ct }. ADMoM amounts to performing the fol- 
lowing iterations (I denotes the iteration index) 

[1.1] Primal update. Given {qj. n (0} an d { v m,nQ)}> update 
primal variables in a coordinate descent fashion; i.e., for n = 
1, . . . , N, update x n as: 



x„(Z + 1) := min £„(x„,Z) 
x„e.4„ 



(19a) 



C n (x n ,l) := £(xi(/), . . . ,x„_i(7),x„,x n+ i(Z), . . . ,x N {l), 

{q J -,n(0}»{«m,n(0}) ( 19b ) 

where £ n (x„, I) is obtained by keeping {x m (l)} m ^ n fixed to 
their values at iteration I. 



[1.2] Dual update. Given the primal variables {x ra (7 + 1)}^ =1 , 
updated multipliers as: 



q,-, n (Z + 1) = q,>(0 + P [*j^n(l + 1) 



-+n,„(Z + l)] 

(19c) 

.(/ + !)] d9d) 



«m,n(' + 1) = Um,n(0 + + 1) _ 

where (3 > is the step-size. 

Once the primal iterates of the neighboring nodes {Xj_> n (1+ 
1)} and {v m {l + 1)} become available at node U n , the dual 
updates d 1 9cb -( fl~9db can be performed locally. As for the 
primal update, the local augmented Lagrangian [cf. (fT8l l 



£ n (x„,Z) = -Un(pn) +C n (P n ) 

c i 



+ 



EA/"„_> 



2 K 



.(0 ~3Cj 



i(OI 



>n,n 1 1 2 



(0 



0' 



< P (0) (20) 

p|raGl p 

can be minimized at node t/„ upon collecting x^j^l) and 
multipliers {q n .i(Z)} from the one-hop neighboring nodes 
Ui G and {w n . P (0} an d ^n,p(0 from nodes C/ p G 

{C^p|C/n S 2p}; that is, from the nodes whose transmissions 
can collide with the ones of C/ n . Roughly speaking, the latter 
quantities pertain to the two-hop neighborhood of node U n and 
are due to the basic properties of the random access strategy. 
If a different medium access protocol such as, e.g., CSMA 
is employed, d!7eb will not be required and the message- 
passing overhead can be further reduced. The ADMoM-based 
distributed algorithm is tabulated as in Algorithm 1, and 
the convergence to the optimal primal arguments {x„(7)} as 
I —> oo is summarized in the following proposition. 

Proposition 1. If there exists a non-zero probability (possibly 
multi-hop) directed path connecting U n to nodes U m G 
Af-^n UA/"„-j. Ul„, for all n, the iterates {x„(Z)} generated by 
Algorithm 1 converge to a globally optimal solution to (P4^' ). 

Proof. Existence of a path connecting U n to nodes U m G 
Af^n U M n ^ U I n , guarantees a regular exchange of local 
primal variables and multipliers among neighboring nodes. 
Under this assumption, convergence of the primal iterates 
{x„(I)} to their optimal values as I — >• oo can be readily 
established using the result in l35l Prop. 4.2]. □ 

Algorithm 1 can also be implemented in an online fashion. 
The equality constraint violation during the initial iterations of 
the algorithm may induce an initial increase of some queues. 
Thus, an online implementation is feasible if nodes can afford 
such a potential increase in the queue length before reaching 
consensus on the local variables. 

A. Handling the interference constraint via primal decompo- 
sition 

Reconsider now problem (P2^), where the interference 
budgets {t^ ax }^! are not partitioned a priori among CR 
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nodes. As primal variables become feasible only when dual de- 
composition algorithms have converged, utilization of network 
parameters obtained from intermediate iterates can possibly 
lead to violation of the interference constraint. To enforce 
strict PU protection during network operation, the primal 
decomposition technique is invoked here; see, e.g., l33ll . 
With this method, resources shared among CR nodes are 
essentially allocated by a master problem. Specifically, at each 
iteration k = 1, 2, . . . of the primal decomposition algorithm, 
problem (P3^' fc )) is solved for given {{t^ a ^(^, k)} Ryn }; then, 
{{ L R a n(t-> k)}n,n} are updated by solving the following prob- 
lem: 

(P5^ fc >) {^(£, k)} = arg m ? 9 ({Q) (21a) 
subject to ig a „ x > 0, (21b) 

2 < V R. (21c) 

To solve (P5^ fc '), the subgradient algorithm can be em- 
ployed EH- Specifically, the subgradient of g{{ k)}) 
with respect to L R la ^(£,k) is given by the optimal La- 
grange multiplier UR, n (k) corresponding to the constraint 
/ lH , n (x„) < L^(£,k) in (P3^) at iteration k J33]. There- 
fore, is updated as 

k + l)= P ls k) + C(fc + l)u R , n (k)} (22) 

where £(•) is the step size, and Vi R {-} denotes projection onto 
the region defined by d21cl i, operation that can be efficiently 
computed as in, e.g. 0361 . At each step k of the primal 
algorithm, CR nodes can employ variables obtained from 
(P3(-W) for network operation, as PU interference protection 
is enforced by updates d22l . 

The projection in (1221 needs to be performed by a "head 
node" in the CR sub-network CR Mr, which is formed by 
nodes that are coupled by constraint 2~2neN /tR,n( x n) — 
L R a nW- P er iteration k, the head node has to collect the 
optimal Lagrange multipliers from the CRs in Mr, and then 
broadcast the updated interference budgets {t^ a ^ c (^, k + l)} n . 
This leads to a semi-distributed algorithm, but the high 
message-passing overhead entailed by centralized solutions is 
nonetheless alleviated. The online algorithm obtained through 
the successive convex approximation and primal decomposi- 
tion is tabulated in Algorithm 2, and its convergence properties 
are summarized next. 

Proposition 2. If there is a cycle connecting nodes U n € Mr, 
V R, the iterates generated by Algorithm 2 converge to a KKT 
solution to (P2). 

Proof. Since the original problem (P2^ ,fc )) is convex, the 
subproblems (P3^ ,fe ^) as well as the master problem (P5^' fc ^) 
are all convex, and thus the globally optimal solution of 
(P2( £ )) is attained via primal decomposition 11331 . Existence 
of a cycle connecting nodes U n E Mr ensures that the 
multipliers {£(&)} can De collected to a cluster head node, 
and that {^^(i, k + 1)} can be subsequently sent back. 
Finally, since (P2^) is optimally solved per iteration I of the 
successive convex approximation, convergence of Algorithm 2 



Algorithm 2 Overall on-line algorithm for (P2) 

Assumption: Path connecting all U„ £ Mr, for all PUs R = 

1, • ■ • , Nr. 

for i = 1, . . . (repeat until convergence) do 

Use solution of (P2<- e - iy ) to compute (l3j-([T4}. If £ = 1, use 
suitable feasible point. 

for k = 0, 1, . . . (repeat until convergence) do 
Receive UR, n (fc) from head node. 
Solve (P3 {e ' k) ) using Algorithm 1. 

Transmit multiplier UR tTl (k) to head node via neighboring 
nodes. 

If head node: update {i%*%(£, k + 1)} via p2). 
Utilize x n (fc) for network operation, 
end for 
end for 



to a KKT point of (P2) is guaranteed (30). □ 

Remark 4 (Fully distributed algorithm). At the expense of 
possibly sacrificing optimality of the resultant exogenous rates, 
coefficients {t/J^} can be set a priori based on the distance 
between CRs and PU R. This may be reasonable especially 
if shadowing can not be estimated ll24ll . In this case, it is not 
necessary to compute the primal decomposition iterates. □ 

V. Packet deliverability in dynamic CR 

ENVIRONMENTS 

Statistics of the SINR may vary during network operation, 
because of the dynamic nature of shadow fading l24l . and the 
variable PU interference levels [cf. (Q])]. CR topology may 
also change with time. Proximity of PUs with intermittent 
activity, or, mobile PU devices may loose link connectivity 
during certain time intervals. The routing problem (P2) must 
be (re-)solved whenever network topology and SINR statistics 
change. Alternatively, it can be implemented online to track 
slow environmental dynamics. Either way, it is necessary 
to establish conditions ensuring that packets are eventually 
delivered to the sink when routes, MAC, and physical layer 
parameters are regularly updated. 

Let s n (r) € {0, 1} be a binary variable taking value 
1 if a packet, after having been randomly routed through 
the network, is placed in U n 's queue at time r, and let 
$„(t) := Pr{s„(r) = 1} denote the probability of such an 
event. Further, collect {i9ti(t)} in the (JV + 1) x 1 vector 
t?(t) := [0i (t), . . . ,0jv+i(t)] t . CR-PU hierarchy may pre- 
vent CR nodes from forwarding packets during certain time 
intervals. Let t n ->j be a binary variable that takes value 1 if 
link U n — > Uj is active, and define Xn->j : = Pr{^n->j = !}■ 
If active, link U„ — > Uj is characterized by a link reliability 
'V-»-.j(' r )' Probabilities {x n ->j} clearly depend on PU activity 
factors and locations, and determine the average connectivity 
of the CR network ifTTI . 

If a packet is in U n 's queue at time r, then U n may 
decide with probability t n -yj (r) to route it through one of the 
available links, where index r emphasizes the time-variability 
of routes. Clearly, if neither node locations nor the PU inter- 
ference or channel conditions change for a certain number of 
time slots, then (t)} and {r„^j(r)} remain invariant. 

The evolution of {^(t - )} can thus be fully characterized 
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by the product probabilities (r) r n ^j(r)}, and the link 

availability factors {% n _^-}. Upon invoking the law of total 
probability, it holds that i? n (r + 1) = Ett 1 Pr i>n( r + 1) = 
l| Si (r) = l,l^ n = l}Pr{ Si (r) = l}Pr{4_>„ = 1} = 
ES *n-tiW' , n->j(T)Xr.-ti'?i(r). Define the (N + 1) x 
(N + 1) packet delivery probability matrix D(t), whose off- 
diagonal entry (i, n) is {t„_>jr„_>iX n _>.i} if t/i is a one-hop 
neighbor of J7 Tl , and otherwise. The diagonal entry (n, n) 
of D(t) represents the probability that a packet remains in 
U n 's queue, which equals 1 - J2i^ n tn-yir n -yiXn->i- Finally, 
since the sink node will not route packets to any other node, 
set the (n,N + l)-th entry of D(t) to D TCj jv+i(t) = 0, and 
Dn+i.n+i( t ) = 1- Matrix D(r) is by construction a column 
stochastic, meaning that D t (t)1at +1 = ljv+i for all r. Then, 
the evolution of {"^(t)} can be expressed in matrix-vector 
form as #(r + 1) = D(t)#(t). 

Using an inductive argument, it is possible to show that the 
(i,n)th entry of the stochastic matrix D(t) := Ilt=i ■^( r ) 
represents the probability that a packet generated at U n reaches 
node Ui in t time slots 1371 Ch. 2]. Therefore, it readily follows 
that a packet is eventually delivered to the sink node U n+i if 
and only if 



lim i?(t) = lim D(i)#(0) 



(23) 



holds for any initial distribution i9(0). A simple condition on 
the CR network topology is provided next in order for (1231 to 
be satisfied. 

Proposition 3. If Ei^nXi-m > 0, V {U n }^ =l and 
En Xn^N+i > 0, a packet stochastically routed according to 
probabilities {t n ^i(r)} over links with reliabilities {r n _>j(r)} 
will be eventually delivered to the destination with probability 
(w.p.) 1. 

Proof. The conditions of Prop. [5] ensure that there exists a 
multi-hop path connecting each node to the destination Un+i 
in the average connectivity graph, where link U n Uj is 
present if Xn-^-j > ifTTl . Let t* be the minimum number of 
time slots such that U n+i can be reached from any node with 
non-zero probability; i.e., t* = mm{t : DM+i, n {t) > V?i = 
1, . . . , N}. Then, the probability that a packet is in Un+i's 
queue at time t* + 1 is given by 

N 

N+l,n(t*)^n(t*) + &N+l(t*) ■ (24) 



Arguing by contradiction, suppose that 
]im.t*-y+oo$N+i(t*) = a < 1; meaning that the packet 
is not delivered to Un+i w.p. 1 — a > 0. Taking the limit on 
both sides of (1241 . one arrives at 



lim 0jv+i(t* + l) 

t*— f +oo 



lim 

t*-S-+oo 



N 



D N+hn {t*)ti n (t*) + t?jv+i(f) 

N 

>min{D N+1 . n {t*)}Y lim i?„(f) + a. (25) 



But since E n =i li m **-s-+oo $n(t*) = 1 — a > and 
min„{_DAr +1 lTi (i*)} > 0, (|25] l can not hold, thus completing 
the proof. □ 
Requiring the existence of an average node-to-destination 
multihop path is tantamount to having a Markov transition 
matrix with a unique absorbing state (the sink node) 
corresponding to the average graph [cf. (l23ll. If a node U„ is 
able to receive packets, but cannot forward them to any other 
node due to a persistent activity of PU nodes in its proximity 
(which violates the condition of Proposition 3), then the 
constraints fi_>„ = for all i G A/l>„ should be added in (PI). 



VI. Numerical results 

Consider the scenario depicted in Fig. [2] where N = 7 
CR nodes cooperate in routing packets to the destination Us- 
Two PU sources also transmit to their intended receivers with 
power 10 dBW. In order to protect the PU system without 
knowing the locations of the PU receivers, 7 points on the 
boundary of the PUs' coverage regions are selected 1241 . The 
PU interference threshold is set to —80 dBW. The path loss 
obeys the model ||x„ — Xj||~' ? , with rj = 3.5. Log-normal 
shadowing is generated with standard deviation 6 dB, and 
to = 1 is used for the small-scale fading (Rayleigh) 0161 . The 
maximum transmit-power of the CR nodes is set to p™ ax = 
dBW, and the noise power is 10~ 8 W. The SINR threshold 
f„ = —10 dB, and the sum of exogenous rates Yln=i P n 
is maximized; that is U n (p n ) = p n and C n {P n ) = 0, for 
all n = 1,...,N. A larger scale network could also be 
considered, but the conclusion that one could draw do not 
depend on the network size. 

Fig. |2a) depicts the optimal routing probabilities 
obtained by solving (P2) with Algorithm 2. At the first 
iteration £ = 1, a feasible starting point is obtained by properly 
modifying the approach of ll38l to the problem at hand, and 
setting the step-size in ( 1221 equal to 1. It can be seen that 
there is a tendency not to route packets through the "southern" 
region of the network; i.e., through nodes that are closer to the 
PU systems. For example, packets generated by U2 are more 
likely to be routed through links U4 — > Uq and Uq — > U7, 
rather than choosing the shortest path JJ% — > U4, — > U5 — > Us- 
Furthermore, node U5 may decide to send packets to Uq rather 
than attempting direct transmission to Us with considerably 
high probability. This is due to the fact that links starting 
from and ending to U4 and [/ 5 are characterized by a higher 
fading- and interference-induced outage probability, as showed 
in Fig.|2jb). In fact, not only PU interference has a detrimental 
effect on the CR SINRs, but also U2, U4, and [7 5 are confined 
to use a lower transmit-power in order to enforce protection of 
the PU receivers. Notice also that [7 2 may decide to transmit 
to U\ instead of U4 with considerably high probability. On the 
other hand, packets generated by U\ and C/3 are routed through 
1/7 with high probability, which in this case coincides also with 
the shortest path. Interestingly, it is necessary to use the primal 
decomposition algorithm only during the first 5-6 iterations 
out of the total 14 (on average) in the successive convex 
approximation algorithm. In fact, the per-CR interference 
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Fig. 2. Test case 1: routing probabilities {t n ^i} (top); and fading-induced Fig. 3. Test case 2 with the second PU transmitter inactive: routing 
outage probabilities (bottom). probabilities {i n _>i} (top); and fading-induced outage probabilities (bottom). 



levels quickly stabilize around steady-state values, with subtle 
variations for I > 6. 

To verify adaptability of the routing probabilities and link 
reliabilities to the states of the PU systems, consider the case 
of Fig. |3ja), where the same CR network operates only with 
PU 1 present. Compared to Fig. EJa), U4 now forwards an 
increased amount of traffic through node Us, As PU 2 is 
inactive, the outage probability of link U4 — > t/5 is lower in 
this case, as confirmed by Fig.|3jb). Furthermore, [/ 5 can raise 
its transmit-power of 10 dB, which significantly decreases 
the outage probability of link U5 — > Us- As a result, almost 
none of the packets (2%) are sent to Uq. Finally, notice that 
CR Uq now splits its traffic evenly between P5 and P7. The 
average exogenous traffic rates, averaged over 20 different 
experiments, are reported in Table [I] It can be seen that p2, 
P4, and P5 increase in this case. This example demonstrates 
the capability of the proposed routing approach to adapt routes 
and transmit-powers to locations of active PUs. 

Fig. |4] corroborates the convergence of Algorithm 1 for (3 = 
0.1 and c € {1, 10}. Specifically, the depicted evolution of 
\t„^j(k)—t n ^j,j(k) \ for nodes U3 and Uq shows that the local 
routing probabilities approximately coincide with those of the 
neighboring nodes after a few iterations. For example, a gap 
smaller than 1% is obtained after 8 iterations. A similar trend 



was observed for the transmit-probabilities, which suggests 
that an online implementation of the algorithm is feasible, and 
queues will be stable after just a few iterations. 



TABLE I 
Exogenous traffic rates. 







u 2 


U 3 


U A 


U 5 


U 6 


U 7 


Test 1 


0.052 


0.023 


0.028 


0.01 


0.011 


0.008 


0.044 


Test 2 


0.057 


0.038 


0.06 


0.05 


0.015 


0.011 


0.051 



VII. Conclusions 

A novel cross-layer optimization framework was intro- 
duced in this paper. Based on channel and interference level 
statistics, and the situational awareness provided by spectrum 
sensing schemes, the novel approach yields optimal routes, 
transmission probabilities, and transmit-powers. The relevant 
optimization problem turned out to be non-convex and hence 
difficult to solve even in a centralized setup. Nevertheless, 
a successive convex approximation was pursued to find a 
KKT solution. Primal decomposition and AD-MoM were 
employed to derive a distributed algorithm, suitable for large 
networks, and amenable to online implementation. As packets 
are randomly routed through the network, their deliverability 
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Fig. 4. Convergence of Algorithm 1 . 



in case of time-varying routing strategies and link reliabilities 
was asserted. Finally, numerical tests verified the ability of the 
proposed scheme to adapt network operation to the propaga- 
tion environment. 



References 

[1] J. L. Wang and J. A. Silvester, "Maximum number of independent paths 

and radio connectivity," IEEE Trans. Wireless Commun., vol. 41, no. 10, 

pp. 1482-1493, Oct. 1993. 
[2] M. Haenggi, "On routing in random Rayleigh fading networks," IEEE 

Trans. Wireless Commun., vol. 4, pp. 1553-1562, Jul. 2005. 
[3] D. D. Couto, D. Aguayo, J. Bicket, and R. Morris, "A high-throughput 

path metric for multi-hop wireless routing," in Proc. Int. ACM Conf. 

Mobile Computing, Networking, San Diego, CA, Sep. 2003, pp. 134— 

156. 

[4] H. Liu, Z.-L. Zhang, J. Srivastava, and V. Firoiu, "Pwave: A multi-source 

multi-sink anycast routing framework for wireless sensor networks," in 

Proc. of Intl. Conf. on Networking, Atlanta, GA, May 2007. 
[5] A. Ribeiro, Z.-Q. Luo, N. Sidiropoulos, and G. B. Giannakis, "Modelling 

and optimization of stochastic routing for wireless multihop networks," 

in Proc. IEEE Int. Conf. on Computer Commun., Anchorage, AK, May 

2007, pp. 1748-1756. 
[6] A. Ephremides, "Energy concerns in wireless networks," IEEE Trans. 

Wireless Commun., vol. 9, no. 4, pp. 48-59, Aug. 2002. 
[7] Q. Zhao and B. M. Sadler, "A survey of dynamic spectrum access," IEEE 

Signal Processing Magazine, vol. 24, no. 3, pp. 79-89, May 2007. 
[8] H. Khalife, S. Ahuja, N. Malouch, and M. M. Krunz, "Probabilistic path 

selection in opportunistic cognitive radio networks," in Proc. of IEEE 

Glob. Telecom. Conf, New Orleans, LO, Dec. 2010. 
[9] C. Xin, L. Ma, and C.-C. Shen, "A path-centric channel assignment 

framework for cognitive radio wireless networks," Mobile Net. Appl, 

vol. 13, no. 5, pp. 463^176, Oct. 2008. 
[10] I. Pefkianakis, S. Wong, and S. Lu, "SAMER: Spectrum aware mesh 

routing in cognitive radio networks," in Proc. of IEEE DySPAN , Chicago, 

IL, Oct. 2008. 

[11] A. Abbagnale and F. Cuomo, "Connectivity-driven routing for cognitive 
radio ad-hoc networks," in Proceedings of IEEE SECON, Boston, MA, 
2010. 

[12] K. R. Chowdhury and I. F. Akyildiz, "CRP: A routing protocol for 

cognitive radio ad hoc networks," IEEE J. Sel. Areas Commun., vol. 29, 

no. 4, pp. 794-802, Apr. 2011. 
[13] S.-J. Kim, E. Dall'Anese, and G. B. Giannakis, "Cooperative spectrum 

sensing for cognitive radios using Kriged Kalman filtering," IEEE J. Sel. 

Topics Sig. Proc, vol. 5, pp. 24—36, Feb. 2011. 
[14] E. Dall'Anese, J. A. Bazerque, and G. B. Giannakis, "Group sparse 

Lasso for cognitive network sensing robust to model uncertainties and 

outliers," Elsevier Physical Communication, Nov. 2011. 
[15] C. Lott and D. Teneketzis, "Stochastic routing in ad-hoc networks," IEEE 

Trans. Auto. Contr., vol. 51, no. 1, pp. 52-70, Jan. 2006. 
[16] G. L. S tuber, Principles of Mobile Communication, 2nd ed. Boston, 

MA: Kluwer Academic Publishers, 2001. 



[17] E. Dall'Anese, S.-J. Kim, G. B. Giannakis, and S. Pupolin, "Power 
control for cognitive radio networks under channel uncertainty," IEEE 
Trans. Wireless Commun., vol. 10, pp. 3541-3551, Dec. 2011. 

[18] L. F. Fenton, "The sum of lognormal probability distributions in scatter 
transmission systems," IRE Trans. Commun. Syst., vol. 8, no. 1, pp. 
57-67, Mar. 1960. 

[19] S. Shin, S. Choi, H. S. Park, and W. H. Kwon, "Packet error rate analysis 
of IEEE 802.15.4 under IEEE 802.1 lb interference," in Proc. of WWIC 
2005, LNCS, Springer, May. 2004, pp. 279-288. 

[20] A. A. Abu-Dayya and N. C. Beaulieu, "Comparison of methods of 
computing correlated lognormal sum distributions and outages for digital 
wireless applications," in Proc. IEEE Veh. Tech. Conf, May 1994, pp. 
175-179. 

[21] R. Rao and A. Ephremides, "On the stability of interacting queues in a 
multi-access system," IEEE Trans. Info. Theory, vol. 34, pp. 918-930, 
Sep. 1988. 

[22] M. J. N. L. Georgiadis and L. Tassiulas, "Resource allocation and cross- 
layer control in wireless networks," Found. Trends in Netw., vol. 1, no. 1, 
pp. 1-144, 2006. 

[23] R. Loynes, "The stability of a queue with non-independent interarrival 

and service times," Mathematical Proc. of the Cambridge Philosophical 

Society, vol. 58, pp. 497-520, 1962. 
[24] E. Dall'Anese, S.-J. Kim, and G. B. Giannakis, "Channel gain map 

tracking via distributed Kriging," IEEE Trans. Veh. TechnoL, vol. 60, 

no. 3, pp. 1205-1211, Mar. 2011. 
[25] R. Zhang, "On peak versus average interference power constraints for 

protecting primary users in cognitive radio networks," IEEE Trans. 

Wireless Commun., vol. 8, no. 4, pp. 2112-2120, Apr. 2009. 
[26] A. G. Marques, N. Gatsis, and G. B. Giannakis, "Optimal cross-layer 

design of wireless fading multi-hop networks," in Cross Layer Designs in 

WLAN Systems, N. Zorba, C. Skianis, and C. Verikoukis, Eds. Leicester, 

UK: Troubador Pub., 2011. 
[27] A. Ribeiro, "Wireless cooperative communications and 

networking," PhD thesis, University of Minnesota. 2006. [Online]: 

http : / /www . seas . upenn . edu/~ar ibeiro/preprints /. 
[28] M. Chiani, D. Dardari, and M. K. Simon, "New exponential bounds 

and approximations for the computation of error probability in fading 

channels," IEEE Trans. Wireless Commun., vol. 2, no. 4, pp. 840-845, 

Jul. 2003. 

[29] N. Ermolova and S.-G. Haggman, "Simplified bounds for the comple- 
mentary error function; application to the performance evaluation of 
signal processing systems," in Proc. of the 1 2th European Signal Proces. 
Conf, Vienna, Austria, Sep. 2004. 

[30] B. R. Marks and G. P. Wright, "A general inner approximation algorithm 
for nonconvex mathematical programs," Oper. Res., vol. 26, no. 4, pp. 
681-683, Jul. -Aug. 1978. 

[31] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge 
University Press, 2004. 

[32] A. Ribeiro, N. Sidiropoulos, and G. B. Giannakis, "Optimal distributed 
stochastic routing algorithms for wireless multihop networks," IEEE 
Trans. Wireless Commun., vol. 7, no. 11, pp. 4261—4272, Nov. 2008. 

[33] D. P. Palomar and M. Chiang, "A tutorial on decomposition methods 
for network utility maximization," IEEE J. Sel. Areas Commun., vol. 24, 
no. 9, pp. 1439-1451, 2006. 

[34] A. Nedic and A. Ozdaglar, "Approximate primal solutions and rate 
analysis for dual subgradient methods," SIAM J. Optim., vol. 19, no. 4, 
pp. 1757-1780, 2009. 

[35] D. P. Bertsekas and J. N. Tsitsiklis, Parallel and Distributed Computa- 
tion: Numerical Methods. Englewood Cliffs, NJ: Prentice-Hall, 1989. 

[36] C. Michelot, "A finite algorithm for finding the projection of a point 
onto the canonical simplex of M n ," J. Optim. Theory Appl, vol. 50, 
no. 1, pp. 195-200, Jul. 1986. 

[37] M. E. J. Newman, Networks: An Introduction. Oxford University Press, 
2010. 

[38] M. Avriel and A. C. Williams, "An extension of geometric programming 
with applications in engineering optimization," Journal of Engineering 
Mathematics, vol. 5, no. 3, pp. 187-194, Jul. 1971. 



IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS (TO APPEAR) 



II 



Emiliano Dall'Anese (S'08, M'll) received the 
^^^^ Laurea Triennale (B.Sc degree) and the Laurea 

/^^^V Specialistica (M.Sc degree) in Telecommunications 

Engineering from the University of Padova, Italy, 
in 2005 and 2007, respectively, and the Ph.D in 
Information Engineering at the Department of In- 
formation Engineering (DEI), University of Padova, 
Italy, in 2011. From January 2009 to September 
2010 he was a visiting scholar at the Department 
of Electrical and Computer Engineering, University 
of Minnesota, USA. He is currently a post-doctoral 
associate at the Department of Electrical and Computer Engineering, Univer- 
sity of Minnesota, USA. 

His research interests lie in the areas of statistical signal processing, 
communication theory, and networking. Current research focuses on wireless 
cognitive radio systems, IP networks, and power distribution networks. 





Georgios B. Giannakis (F'97) received his Diploma 
in Electrical Engr. from the Ntl. Tech. Univ. of 
Athens, Greece, 1981. From 1982 to 1986 he was 
with the Univ. of Southern California (USC), where 
he received his MSc. in Electrical Engineering, 
1983, MSc. in Mathematics, 1986, and Ph.D. in 
Electrical Engr., 1986. Since 1999 he has been a 
professor with the Univ. of Minnesota, where he now 
holds an ADC Chair in Wireless Telecommunica- 
tions in the ECE Department, and serves as director 
of the Digital Technology Center. 
His general interests span the areas of communications, networking and 
statistical signal processing - subjects on which he has published more 
than 300 journal papers, 500 conference papers, 20 book chapters, two 
edited books and two research monographs. Current research focuses on 
compressive sensing, cognitive radios, network coding, cross-layer designs, 
wireless sensors, social and power grid networks. He is the (co-) inventor of 
twenty patents issued, and the (co-) recipient of eight paper awards from the 
IEEE Signal Processing (SP) and Communications Societies, including the G. 
Marconi Prize Paper Award in Wireless Communications. He also received 
Technical Achievement Awards from the SP Society (2000), from EURASIP 
(2005), a Young Faculty Teaching Award, and the G. W. Taylor Award for 
Distinguished Research from the University of Minnesota. He is a Fellow of 
EURASIP, and has served the IEEE in a number of posts, including that of 
a Distinguished Lecturer for the IEEE-SP Society. 



